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Abstract 

A set of action-angle variables for a soliton cellular automaton is obtained. 
It is identified with the rigged configuration, a well-known object in Bethe 
ansatz. Regarding it as the set of scattering data an inverse scattering method 
to solve initial value problems of this automaton is presented. By considering 
partition functions for this system a new interpretation of a fermionic character 
formula is obtained. 



1 Introduction 



X 



1.1 Action-angle variables for the box-ball system 

The box-ball system is a well-known example of soliton cellular automata [1, 2]. Here 
we show its simplest version. 
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---0 12345678 9 101112 131415 16171819 20 2122 23--- 

The arrows show the time evolution process of the automaton whose precise rule 
will be given in the main text (See Section 2, also for the convention of the position 



coordinate). This system and its generalization to accommodate several kinds of 
balls are known as integrable systems by a certain limiting procedure on discrete 
KdV/Toda equations [3, 4] or by explicit constructions of their conserved quantities 
[5, 6]. Recently Kuniba, Okado, Yamada and the author presented a conjecture [7] 
that a set of action-angle variables for this system was given by the rigged configura- 
tion, an object known in Bethe ansatz [8]. In the simplest case a rigged configuration 
is a Young diagram with a set of integers rigged to its rows. For instance according 
to the above example the associated rigged configuration evolves as 
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The rows of the Young diagrams are regarded as the action variables and the numbers 
(riggings) associated with each row the angle variables. As one sees the former are 
invariant under the time evolution while the latter increase by the lengths of the 
associated rows. In other words the time evolution is linear in these variables. 

In this paper we present a new method to construct a rigged configuration from 
a given state of the box-ball system. Here the procedure to construct a Young 
diagram is the same as that in [6, 9] but our method can extract additional data, 
the riggings. Having the additional data our method provides a bijection between 
the set of automaton states and the set of rigged configurations. To borrow a name 
from the theory of nonlinear evolution equations [10] we call this bijection the inverse 
scattering transform for the box-ball system, where we regard the rigged configuration 
as a set of scattering data. Using the technique of the inverse scattering transform 
we give a proof of the assertion in [7] that the rigged configuration is indeed playing 
the role of a set of action-angle variables. 



1.2 The box-ball system and the fermionic formulas 

We observe that a certain crystallization of solvable lattice models can produce soli- 
ton cellular automata [11, 12, 13], while there is a connection between solvable lattice 
models and Bethe ansatz in a character level [14, 15]. Now the inverse scattering 
transform for the box-ball system is in some sense providing a direct connection 
between soliton cellular automata and Bethe ansatz. The second purpose of this 
paper is to complete a study on this connection in the character level which was also 
proposed in [7]. 



Let us consider a particular case where the character formula comes from two 
different expressions for Kostka polynomials [16]. One is from the Bethe ansatz [17] 
and yields the fermionic character formula which is made of sums of products of 
g-binomial coefficients. The other is from the solvable lattice models in statistical 
mechanics [18] and yields the one-dimensional configuration sums over paths (tensor 
products of crystals). To be more specific we consider the Kostka polynomial K^y{q) 
with /i = {L~s^s) and v = (1^) where L and s are two integers obeying L/2 > s > 0. 
Then the two expressions for the Kostka polynomial lead to the identity [14] 



K^^^iq) = Y^qm-LT.>,^.W+'-^^ JJ 
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bePL,. 



where A h s means that A is a partition of s. Here Pi{X),mi{X) and 0(A) will be 
defined in the main text (equations (11), (12)), and 6'(true) = l,^(false) = 0; the set 
of highest weight paths is defined as 



P- 



L,s 



(fel,...,&L) 



h e {0, 1}, ^bj < - for any i, ^ 6^- = s [> . 



There is a decomposition Pl^s = LJai-s-Pl(A) such that if b G -Pl(A) then the relation 

■13-- 



^ • ^ Q{hj < 6j+i) = J2i>i''^iW holds. Therefore equation (1) leads to 
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Furthermore the decomposition makes each term of the sum XIai-s vanish separately. 
Thus there is actually a set of refined identities behind (1) with respect to A. In the 
context of Bethe ansatz this A has a clear meaning; it is a label for the associated 
eigenstate of the Heisenberg magnet. However in the context of solvable lattice 
models its meaning is so far unclear. By identifying the space of paths for a lattice 
model (the six-vertex model in f err o- magnetic regime) with the space of states of 
the box-ball system, we shall find that the path associated with A = (I™'i2™'2 _ _ /) 
has Tjii "solitons" of length i for any i. This interpretation is based on an expression 
for partition functions for the box-ball system with specified soliton contents, which 
turns out to be the above mentioned refined identities of the fermionic formula. This 



expression for partition functions and its variant for including non-highest weight 
paths were proposed in [7]. We derive these identities (Theorem 34) as a resuh of 
the inverse scattering transform for the box-ball system. 

1.3 Outline of the paper 

In Section 2 the box-ball system is introduced and its updating rule is explained in 
terms of arcs. In Section 3 we define a mapping S which sends any state of the 
box-ball system to a rigged configuration (set of scattering data). The inverse map 
of S (inverse scattering transform) is defined in Section 4. The notion of wave tails 
and wave fronts of automaton states is also introduced here. In Section 5 we give 
a proof that the time evolution is linearized in the scattering data , and present an 
inverse scattering method for this system. 

In Section 6 we derive formulas for the right-most wave front and the left-most 
wave tail of the automaton state associated with a given rigged configuration. In 
Section 7 automaton states in a finite interval are studied and the notion of highest 
weight states is introduced. By taking summations over these states we define par- 
tition functions for the box-ball system in Section 8 and establish the conjectured 
fermionic formulas proposed in [7]. In Section 9 we prove that our inverse scattering 
transform is equivalent to the s/(2) case of the bijection in [8]. 

In Appendix A we give a proof and an explanation of two formulas for one- 
dimensional configuration sums expressed by g-binomial coefficients. 

2 The box-ball system 

We consider a one-dimensional array of infinite number of boxes that extends towards 
both directions. As a position coordinate we put successive integers to the walls 
between the boxes rather than to the boxes themselves. 

Any box is either an empty box or a filled box. The latter means that there is a 
ball within the box. We assume that there are at most finite number of balls in the 
system. The empty box is denoted by FH and the filled box is by % . Clearly there 



are four types of configurations of adjacent boxes, , #, # and 



We adopt the updating rule of the box-ball system in [9] : 



1. For every % connect its two boxes with an arc. 



2. Ignore those boxes connected with the arcs and regard the other boxes as if 
they were successively adjoining. 



3. Repeat steps 1 and 2 as many times as possible. 

4. For every pair of connected boxes interchange % and | | . 
Example 1. We have 





■■0123456789 10 11 0123456789 10 ll^^^ 

Definition 2 (advanced/retarded arcs). The arcs that one will obtain by ap- 
plying the items 1 (resp. 1 but # was replaced by #), 2, and 3 in the above 
procedure to any automaton state are called the advanced arcs (resp. retarded arcs ) 
of that state. 

By definition we have that 

Lemma 3. Let p and p' be two automaton states. Suppose the set of advanced arcs 
of p is equal to the set of retarded arcs ofp'. Then p' is the state that one will obtain 
from p by updating the system once. 

For later use we introduce the notion of 

Definition 4 (depth). If an arc has no arc within it, its depth is set to be one. If 
an arc has arc(s) within it, its depth is given by 1 + maxjdepths of inside arc(s)}. 

For instance there are three depth 1, one depth 2, and one depth 3 arcs in each 
figure of Example 1. 



3 Scattering data and soliton contents 



Given a set of integers {ai}i<i<2n "we introduce a matrix of the form 

M = 



ai as 
a2 c^A 



«2n-l 
a2n 



(2) 



Let Tin (resp. j-tn) be the set of all matrices of this form subject to the condition 
tti < tt2 < ■■■ < tt2n (resp. «! < ^2 < ■ ■ ■ < tt2n)- We define Ti := UJ^qTY^i 

o o o 

(resp. Ti '■= LJ^Q7i„). Here the set Tio = Tio consists of a formal two-row matrix 

o 

with no element. We associate any matrix M G Tin with an automaton state that 
has balls between the walls a2k-i and a2k for 1 < A; < n. We shall occasionally 
identity such a matrix M with the associated automaton state itself. 



Example 5. We have the identification 

14 9 



^ M = . 

0123456789 10 11- ■■ ^^6 10 

We show a procedure that associates a rigged configuration, a Young diagram 

o o 

with a set of integers, for any M E Ti. Suppose M G Tin and let rii = n and 
Ml = M. Given any non-negative integer i and a matrix of the form 

we define a set of integers /i^, r;^, a^, . . . , a^. (tlie upper indices are not exponents) and 

o 

a matrix Mj+i G 7ini+i- Here n^+i = n^ - t;^. We set hi = mmi<k<2n,~iioil+i - al) 
and Pi = Oil — ^^i- From the set of integers {Pl}i<k<2n^ we remove pairs of identi- 
cal integers repeatedly until there is no such pair. Let Vi = # (removed pairs). We 
arrange the removed integers (after their multiplicity was divided by two) in increas- 
ing order and call them a\, . . . ,al.; we call the remaining integers a\^^, . . . , Oi'^2nl+i i^ 
increasing order. Then let Mj+i be a matrix of the form (3) but with all the i's were 
replaced by i + I's. By repeating this we shall obtain a matrix M^+i which has no 
element after a finite number (/) of steps. Then let A be the Young diagram in Figure 
1 and a = LJi<j</{a*}i<j<t,. the set of integers obtained through this procedure. 
We denote TZC (resp. 3^) the set of all rigged configurations (resp. Young dia- 

o o 

grams). Let ^ : 7i ^ "R-C (resp. ^ : 7i ^ 3^) be the mapping that sends any matrix 

o 

M G 7i to a rigged configuration ^(M) = (A, a) G TZC (resp. a Young diagram 
^{M) = A G 3^) by the above procedure. For a reason that will be clear afterwards 
we call them 

o 

Definition 6 (scattering data, soliton content). For any M G 7i (or the automa- 
ton state associated with M) the rigged configuration S{M) is called its scattering 
data, and the Young diagram ^(M) is called its soliton content. 

o 

The mapping ^ becomes a bijection between 7i and 7ZC. Its inverse map will be 
discussed in the next section. 

Example 7. We consider the matrix M in Example 5. Then we have hi = min{3 — 
1, 4 - 3, 6 - 4, 9 - 6, 10 - 9} = 1. Since 

14 9\_/13 5\ . _ ( 1 4 
36 10 J \2A6J^~\12 4 
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Figure 1: The Young diagram. 



we obtain Vi = 2, a\ = I, a^ = A, and M2 = I <-> ) • Clearly /i2 = 2 and 



since 



X 2 



we obtain V2 = 1 and af = —2. Therefore we have 



1 4 9 
3 6 10 



4 
1 



We usually depict a rigged configuration by putting the riggings along the associ- 
ated vertical edges of the Young diagram (See Example 7). This convention matches 
to the following step by step description of the mapping S. Suppose we are on a 
square lattice made of vertices connected by unit length bonds, and have a matrix 

o 

M G 7i. We start at an arbitrary vertex of the lattice. Every step we proceed right- 
ward or upward by one unit length and replace our matrix M by the following rule. 

o 

First suppose M ^ Tio- 



1. If M G 7^ we proceed rightward and replace M by M — (24''^ that is in H. We 
place a horizontal segment of unit length on the way. 



2. \i M E Ti but M ^ T-C there are pairs of same integers in M. Find a pair of 
smallest same integers, one of which is in the first row and the other in the 
second row. We proceed upward and remove the pair from M. We place a 
vertical segment of unit length on the way and put the removed integer on its 
right-hand side. 

We repeat this procedure until our matrix M has no more entry. Then we stop there. 
By connecting the segments we obtain a path of lower-left to upper-right direction. 
By attaching an upper left corner to it (to complete a Young diagram) we obtain the 
rigged configuration. 

4 Inverse scattering transform 

In this section we study the inverse of the mapping S. For any integer x we define 
a map fij; : H — > 7i by the following rule. Given M G Tin we label its elements as 
(2). Let i be the largest integer obeying the condition ai < x; if ai > a; we let i = 0. 
For an odd i(= 2/ — 1) we set 

fix(M) : = 

and for an even i{= 21) 

Q (nj) .= Z' '^1 ■ ■ ■ '^21-1 X a2l+l ■ ■ ■ 0L2n-\ 
\«2 ■■■ 0121 X a2l+2 ■■■ Oi2n 

We extend our definition to any set of integers X. Given M E Ti and X = 
{xi, . . . ,Xp} we set 

o 

In addition we define the following maps ^i : Ti ^^ Ti (i = 0, 1). Given M E Tir. 
we set 

For later analyses it is important to notice that the automaton state associated with 
$1 o Qx{M) (resp. $o ° ^x{M)) is made by the following algorithm. We begin with 
the automaton state associated with M. 



ttl •• 


■ ■ tt2«-l X 


' ■ 0L2n-\ 


a2 ■ ■ 


- X a2i ■■ 


■ ■ «2n 



1. At the wall positions Xi, . . . , Xp and ai, 0:3, . . . , Q;2n-i (resp. 0:2, ^4, • • • , C(2n), 
put marks, say V. Here their multiplicity should be taken into account. 

2. Split the array of boxes at the left-most wall with the mark(s). If it has 
p{> 1) mark(s), translate every box and ball on the right-hand side of the wall 
rightwards by the width of 2p boxes. Delete the mark(s) at the wall. 

3. Repeat item 2 as many times as possible. 



4. Fill 



s (resp. 



3) into the gaps between the arrays of boxes made by 



the above procedures. 

Remark 8. This algorithm also applies to $0 and $1 themselves, because we can 
regard them as $0 ° ^x and $1 o Qx with X = 0. 

Example 9. Let M be the matrix in Example 5 and let X = {1,5,8}. Then we 
have 



$iol]^(M) 



1 1 A 5 8 9 

1 3 5 Q 8 10 

2 4 9 12 17 20 

3 7 11 14 18 22 



13 5 7 9 11 
2 4 6 8 10 12 



and 



$oofix(M) 



1 3 8 11 16 19 

2 6 10 13 17 21 



The algorithm to make the automaton state associated with $1 o Qx{M) is depicted 

as 

V 

V V V V 



^ 
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0123456789 10 1112 13 1415 16 1718 19 20 2122- 



In contrast that for $o ° ^x{M) is depicted as 

V V V V V V 



0123456789 1011- 

i 
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0123456789 10 1112 13 14 15 16 1718 19 20 2122- 



Definition 10 (wave tails/fronts). In the state of the box-ball system associated 

o 

with M G Ti the wall positions specified by the first row of M are called wave tails. 
The wall positions specified by the second row of M are called wave fronts. 

Remark 11. For any set of integers X the mapping $i o Qx (resp. $0 ° ^x) embeds 



the box pair (resp. % ) at every wave tail (resp. wave front) of M, and 
at every wall position specified by X. This implies that the set of retarded arcs 
of $1 o flx{M) (resp. advanced arcs of $0 ° ^x{M)) consists of — 1) depth > 2 
arcs made out of the retarded arcs (resp. advanced arcs) of M, being stretched but 
their topology unchanged, and — 2) depth 1 arcs introduced along with the newly 
embedded box pairs. See Figure 2 which corresponds to the latter case in Example 
9. 

Let (A, a) be a rigged configuration where A is the Young diagram in Figure 1 and 
a = \-\i<i<i{a'j}i<j<v. a set of integers. For any i we denote {a]}i<j<Vi by a* (Again, 

o 

the upper indices are not exponents). Let L/+i(A,a) = M;+i(A,a) G Ho be the 
two- row matrix that has no element. For any i (1 < i < /) we define the matrices 

Li{\, a) , Mi{\, a) e Hn, (where n, = Y.k=i^k) recursively as 

Li{\,a) := ($0)^' ofi^.(Li+i(A,a)), 

M,(A,a):=(<l>i)^>of]„,(M,+i(A,a)). 

We write L{X,a) := Li(A,a) and M(A, a) := Mi(A,a). Then we can regard L and 

o 

M as mappings from TZC to 7i- It is easy to see that the map M is the inverse of 
the map S in Section 3. 
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12 3 4 5 6 7 



9 10 11 12 13 14 15 16 17 18 19 20 21 22- 



Figure 2: The embedding of box pairs and the advanced arcs. 

Proposition 12. S{M{\,a)) = (A, a). 

We call this map 
Definition 13 (inverse scattering transform). For any (A, a) G TZC the mapping 

o 

(or its image) (A, a) i-^ M{X, a) G 7i is called its inverse scattering transform. 

The inverse scattering transform and the mapping S (which may be called direct 

o 

scattering transform) form a bijection between TZC and H- We will also call this 
bijection itself the inverse scattering transform for the box-ball system. 



5 A proof of the linearized time evolution of the 
scattering data 

The generic appearance of the automaton states in the real space changes in not 
a simple way under the time evolution. According to the time evolution the set 
of associated scattering data also evolves. In this section we show that the time 
evolution of the scattering data is rather simple; it turns out to be linear. 

Recall the two-row matrices Mj(A, a) and Lj(A, a) in Section 4. For any i we have 

11 



Lemma 14. The first row of Mi{X,a) is equal to the second row of Li{X,a). 

Proof. We give a proof by descending induction on i. For i = I + 1 the claim holds 
since the matrices have no entry. Let {Pi, . . . , Pp) be the first row of Mj_|_i(A,a) 
and suppose it is equal to the second row of Li^i{X,a). Then the first row of 
Oai(Mj_|_i(A, a)) is equal to the second row of r2ai(Lj+i(A, a)) since both of them are 
made out of {a\, . . . , a^,, /3i, . . . , Pp} by rearranging them in increasing order. The 
claim of the lemma follows immediately because of the definitions of the mappings 
$0 and $1. D 

For any i we have 

Lemma 15. The set of advanced arcs of Li{X,a) is equal to the set of retarded arcs 
ofM,{X,a). 

Proof. We give a proof by descending induction on i. For i = I -\- 1 the claim holds 
since the matrices have no entry. Suppose that the claim is true if i was replaced by 
i + 1. First we will show the following assertion: 

• The set of advanced arcs of $0 ° ^a^{Li+i{X, a)) is equal to the set of retarded 
arcs of $1 o r2a»(Mj+i(A, a)). 



Here the <l>o°^a» embeds # 3 to the wave fronts of Lj+i(A, a) and the wall positions 



specified by the set a*; the $1 o f2„i embeds #3 to the wave tails of Mi^i{X,a) 
and the wall positions specified by the same a*. But the wave fronts of Lj+i(A,a) 
are the wave tails of Mj+i(A,a) by Lemma 14. Therefore every embedding point 
coincides exactly between both cases. This together with the assumption of the 
induction concludes that the above assertion is true (See Remark 11). By repeating 
this argument we can obtain the claim of the lemma for i (See Remark 8). The proof 
is completed. D 

By Lemma 3 it leads to 

Corollary 16. // the state of the box-ball system associated with L{X, a) is updated 
once, it is the state of the system associated with M{X,a). 

We give an example for Lemmas 14 and 15. 
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Example 17. Recall the second rigged configuration of the example in Section 1. 



(A, a) 



The A has the shape of hi = h2 = Ijh^ = 3 and vi = 2,V2 = v^ = 1. The set of 
riggings is a = a^ L\ a^ L\ a^ with a^ = {7, 9}, a^ = {2}, a^ = {0}. Thus we have (by 
denoting Li{\,a) by Lj) 

$0 o n,s (L4) = (^^ , <^l o fi,3 (L4) = ("2) , L3 = <^lo n,s (L4) = Q , 
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L2 = $oofia2(^3) 



0^\ (0 2\_/0 4 
^ Syl^Vl 3y'~V3 6 

4 7 9\ /0246 
3 6 75J^ll35 7 



6 11 15 
4 9 12 16 



Their advanced arcs are given as follows. 
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On the other hand we have (by denoting Mj(A, a) by Mj) 

$l0f],3(M4)= Q), $?ofi,3(M4)= QY M3 = $?ofi,3(M4)= fg', 



M2 = $1 oQ^2{M3 



Ml = $1 ol]„i(M2 



^ 3 
^ 6 



3 6 7 5 



1 3 

2 4 



3 6 

4 10 



13 5 7 



4 75 10yV246 
Their retarded arcs are given as follows. 



4 9 12 16 
6 11 15 18 
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0123456789 101112131415 16 171819 20 2122 23- 
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Let a be the rigging made out of a by replacing a* by a* := a* — X]l=i ^fe- Let a* 
denote the set {a]}i<j<vi- We set Ti := ^q""*^ o $;^. By the same symbol we write as 
Ti(X) = {xi + 1, . . . , Xp + 1} for any set of integers X = {xi, . . . , Xp}. For instance 
we have a* = T-^ ^~^ '' '(a*). It is easy to see that 

Lemma 18. For any M ETi and any set of integers X the relation VLtj^{^x){Ti{M)) = 
TioVLx{M) holds. 
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Then for any i we have 

Lemma 19. The following relation holds: Li{X,a) = T^^ '"^(Mj(A, a)). 

Proof. We give a proof by descending induction on i. For i = I + 1 the claim 
holds because both sides are equal to a matrix with no entry. Suppose Lj+i(A,a) = 
T-^^^'"^ '(Mj+i(A, a)). By repeated use of Lemma 18 we have J7ai(Lj_|_i(A, a)) = 
T^+-+^^of]„.(Mi+i(A,a)). Then 

Li{\,a) = <l>o'ofi„.(L,+i(A,a)) 

= ($0-^ o $,)^i+-+'*-i o $^ o fi,.(M,+i(A, a)) 



t; 



i/»lH h/l-i- 



^(M,(A,a)). 



D 



By setting i = 1 in this lemma we have 
Corollary 20. The following relation holds: L{X,a) = M{X,a). 

We give an example for Lemma 19. 
Example 21. Recall the first rigged configuration of the example in Section 1. 



(A, a) 


















8 
6 





-5 



The set of riggings is a = a} L\ a^ L\ a'^ with a^ = {6, 8}, a^ = {0},(£' = {—5}. We 
have (by denoting Mj(A,a) by Mj) 



Ml = $ior2„i(M2 



$2ofi,3(M4) 



-2 0\ / 1 3 
1 j ^ I 2 4 



-1 



, M3 = $?ofi,3(M4) 



-2 
1 



-13^ 
2 5^ 



+ 



-1 3 
2 5 

13 5 7 



6 11 15 
4 9 12 16 



2 4 6^ 
Compare the M3, M2, Mi in this example with the L3, L2, Li in Example 17. 
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Now we present the main result of this paper. By combining Corollaries 16 with 
20 we have 

Theorem 22. // the state of the box-ball system associated with M{X,a) is updated 
once, it is the state of the system associated with M(A,a). 

By this theorem we can regard the rows of A as the action variables for the box- 
ball system, and the riggings as the associated angle variables. Initial value problems 
of the box-ball system can now be solved by an inverse scattering method. Let T{a) 

o 

be the rigging made out of a by replacing a* by a* + X]fc=i ^fc- Let T{M) G 7i be 

o 

the matrix for the automaton state which is made out of the state for M G 7i by 
updating once. Then for any integer A^ we have 



Direct Scattering 



Inverse Scattering 



Automaton State 

Men 

Time Evolution 

Automaton State 
T^{M) G n 

This diagram shows that the time evolution of any automaton state is given by a 
composition of the direct scattering, the linearized time evolution of the scattering 
data, and the inverse scattering. 

To close the section we recall the notion of soliton content in Definition 6. We 
can write the soliton content A as a partition 



Scattering Data 


(A,a)G7^C 




Linearized 


Scattering Data 


(A,T^(a))G7^C 



A= (l'"i2'"^../i™''). 



(4) 



Here mi G Z>o is the multiplicity of i in A. We assume rrih 7^ 0; if we take A as 
in Figure 1 then h = Ylk=i^k, namely the width of the Young diagram. We can 
say that there are rrih solitons of length h, rrih-i solitons of length h — 1 , ... , and 
mi solitons of length 1 in the automaton state associated with A. The numbers of 
successive balls in the real space can change. However the solitons (in the space of 
scattering data) keep their independence. To see the soliton content in the real space 
one generally has to apply the updating process many times. 



16 



6 Upper and lower bounds for the inverse scatter- 
ing transform 

Before describing the aim of this section, here we give a representation for rigged 
configurations that is shghtly different from the previous one. Let A be the Young 
diagram of the form (4). Denote the set of riggings for i by J* = {Jj}i<j<mi, J] ^ ^• 
In other words J* is associated with the set of A's rows of length i. Then if ttij = we 
have J* = 0. Let J = Ui<j</i J*. Now the pair (A, J) specifies a rigged configuration. 

Remark 23. For reasons of simpler descriptions of some statements we shall assume 
Jl<Ji<---<Jln, ^01 any i. 

In this section we study the upper and lower bounds for the inverse scattering 
transform. The result is used in the following sections to study the relation between 
partition functions for the box-ball system and fermionic formulas. It is also used 
to prove that our inverse scattering transform is equivalent to a bijection in [8] 
(Section 9). 

For any set of integers X we denote $10 Vtx by IIx, where $1 and Vtx are the 
maps introduced in Section 4. We set 

M,(A, J) = II ja o Hj.+i o . . . o nj.(M;,+i), (5) 

M(A,J) = Mi(A,J), (6) 

o 

where M/j+i G Tio is the matrix with no entry. We can regard M as a mapping from 

o 

IZC to 7Y- It should be observed that this M is the same as the M in Section 5 
although the descriptions of rigged configurations are different. By abuse of notation 
we shall denote S""-^ either the M or the M. 

We try to find expressions for the lower right and the upper left elements of the 
matrices Ma(A, J) in (5) because they determine the upper and lower bounds for the 
inverse scattering transform. Given any i with ttij > let c*'" = X]fc=a (1^111(^5 i) — a + 1) ruk 
and J^'"' = Jj + 2d'"- for 1 < a < ? + 1. It is easy to see that 

ja-l,a ^ ja-1 (7) 

and 

c'^^+^ = d^^-{ma + --- + mh). (8) 

For 1 < a < /i we have 



17 



Lemma 24. 

1. The lower right element of Qja[Ma+i{X, J)) is m.aXa<i<h,mi>o{J^~^^} ■ 

2. The lower right element of Ma{X, J) is m.a:Xa<i<h,m^>o{Jni^}- 

Proof. We give a proof by descending induction on a. For a = h the item 1 holds 
because both sides are equal to J^^. Assume that the item 1 holds for a. Since 
MaiX, J) = $1 o nja(Ma+i{\, J)) we have 

L. R. E. of M,(A, J) = L. R. E. of fija(M,+i(A, J)) + 2(m, + ■ ■ ■ + tjih) 

a<i<h,mi>0 
a<i<h,mi>0 

(L. R. E. is for lower right element.) Here the first equality is by the definition of $i, 
the second is by the assumption of induction, and the third is from (8). Therefore 
the item 2 holds for a. Then we have 

L. R. E. omj.-.(Mai\,J)) = max{J^-^\, max {JX}} 

a<i<h,mi>0 ' 

= max {Jm-}- 

a-l<i<h,mi>0 '■ 

Here we used (7). Therefore the item 1 holds for a — 1. The proof is completed. D 

For 1 < a < h we have 
Lemma 25. 

1. The upper left element of Qja(Ma+i{\, J)) is rmna<i<h,m,>o{Ji + i — a}- 

2. The upper left element of Ma(A, J) is rama<i<h,mi>o{Ji + i — a + I}. 

Proof. We give a proof by descending induction on a. For a = h the item 1 holds 
because both sides are equal to J^^. Assume that the item 1 holds for a. By the 
definition of $1 the item 2 for a follows. Then we have 

U. L. E. offija-i(M„(A, J)) = min{Ji"-\ min {Jl + i-a + l}} 

a<i<h,mi>0 

== min {Jl+i - {a- 1)}. 

a—l<i<h,mi>0 

(U. L. E. is for upper left element.) Therefore the item 1 holds for a — 1. The proof 
is completed. D 



Then we have that 

Theorem 26. Let (A, J) be a rigged configuration of the form given at the beginning 
of this section. For the automaton state associated with S~^{X,J) the following 
statements hold. 

1. The number of balls is |A| := X]/i>i>i ^^^j- 

2. The left-most wave tail of the state is mini<j</j„j.>o{<^i+^}; '^'^^ ^^^ right-most 
wave front is maxi<j</,,„^>o{J;,^ + 2X1^=1 min(i,j)mj}. 

Proof. Item 1 follows from the fact that each of the Ilja in (5) is embedding as many 
^s as the height of the a-th column of the Young diagram A, which is equal to 



nia + ■ ■ ■ + nih- Item 2 follows from Lemmas 24 and 25. D 

7 Automaton states in a finite interval 

o 

In the set of all automaton states 7i we define a set of all the states, and a set of all 
the highest weight states in a finite interval with a specified soliton content. Then 
we determine the associated subsets in the set of all rigged configurations TZC. The 
results are used in the next section. 

To begin with we introduce the notion of 

Definition 27 (highest weight state). A state of the box-ball system is called a 
highest weight state if all the balls are on the right-hand side of the wall 0, and for 
any l the number of empty boxes between the walls and d, is not less than that of 
filled boxes there. 

Consider an array of L boxes. We put the numbers 0, 1, . . . , L to its walls. Recall 
the matrix notation in Section 3. Let 7in{L, A) be the set of all matrices of the form 
(2) obeying the conditions < ai,a2n < L and ^(M) = A, and H^{L,\) be its 
subset with additional conditions 2 J2k=i '^2fc + tt2j < 2 J2k=i "^afe-i for I < j < n. In 
the latter case the set of additional conditions is equivalent to that for the highest 
weight states in Definition 27. We set 

n{L,x) = Un>onn{L,x), (9) 

7^+(L,A) = U„>o7^+(L,A). (10) 

They are sets of automaton states in the finite interval between walls and L, and 
with the specified soliton content A. 

19 



Recall the rigged configuration (A, J) of the form given at the beginning of Sec- 
tion 6. Fix a positive integer L. As functions of the variable A we write 

Pi(A) : =L-2^min(i,j)mj, (11) 

mi{\) : = rrii, 0(A) := ^ m.in{i, j)mimy (12) 

The latter will be used in the next section (Lemma 33). We set 

7^C(L, A) = {(/i, J) e nC\ii = A, -i < Jj < ■ • ■ < J^^ < pi{\) for any i}, (13) 
7^C+(L, A) = {(^, J) e nC\i2 = A, < Jj < ■ ■ ■ < J^^ < pi{X) for any i}. (14) 

By definition we have n^{L,X) C H(L, A) and nC+{L,X) C nC{L,X). From 
Theorem 26 it is clear that 

Lemma 28. The inverse scattering transform is a bijection between H{L, A) and 

nc{L,x). 

Moreover we have that 
Lemma 29. The inverse scattering transform is a bijection between 7i^{L,X) and 

nc+{L,x). 

Proof. Recall that the map S~^ is so defined as to make automaton states by em- 



bedding 3 recursively at the wave tails, and at the wall positions specified by 
the riggings. We assume J^, < Pi{X) for any i; if otherwise the right-most wave 
front of S~^{X,J) would exceed L. Suppose (A, J) G 71C~^{L, X). By induction we 
prove the automaton state associated with Ma{X, J) in (5) is a highest weight state. 
li a = h + 1 the claim holds by definition. Suppose the automaton state associ- 
ated with Ma+i{X, J) is a highest weight state. By the embeddings we are shifting 
the divided arrays of boxes to the right (See Example 9). Therefore if a state is a 



highest weight state, then the state obtained from it by embedding a at any 
wave tail, or at any non-negative wall position is a highest weight state. Thus by 
the assumption of induction the automaton state associated with Ma{X, J) is also 
a highest weight state. By taking a = 1 we have S^^{X,J) G 1-C^{L,X). Thus 

s-\nc+iL,x))cn+iL,x). 

Conversely suppose (A, J) ^ 1ZC^{L, A). Then there is a negative rigging J{{< 0) 



for some i. By embedding a % sX J{ we obtain a state that is not a highest weight 
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state. One cannot obtain a highest weight state from a non-highest weight state by 



embedding a % at any position. Therefore we have ^ ^(A, J) ^ 7i^{L, A). Thus 



s{n+{L,x))cnc+{L,x). n 

8 Partition functions with a specified soliton con- 
tent 

In this section we study partition functions for the box-ball system. Recall the 
quantity rrii in (4). In the context of the box-ball system the m, was the number 
of solitons of length i. By considering the partition function taken over the highest 
weight states, we shall find that the nii coincides with a well-known quantity in 
the context of Bethe ansatz, the number of strings of length i. We also consider a 
partition function taken over not necessarily highest weight states. 

First we introduce the notion of energy for states of the box-ball system. Among 

the four types of configurations of adjacent boxes, let > # # > # do not cost 



the energy, and let % cost the energy of its center wall position. 



Definition 30 (energy of an automaton state). For any matrix M G 7i of the 
form (2) (or for the automaton state associated with M), the quantity EcrmiM) : = 
^k>i '^2fe-i is called its energy. 

This energy is not a conserved quantity of the box-ball system. It has its origin 
in the corner transfer matrix (CTM) analyses of the solvable lattice models [19]. In 
particular the energy in Definition 30 is for the six-vertex model [20] but the sign 
was reversed according to the change of the model to its ferromagnetic regime. 

Next we introduce the notion of 

Definition 31 (energy of a rigged configuration). For any rigged configura- 
tion (A, J) of the form at the beginning of Section 6 the quantity E^c{^, J) '■ = 
T.h>i,j>i vcim{i, j)mimj + X^jLi EJ^i Jj is called its energy. 

We observe that the inverse scattering transform is an energy-preserving bijection 
in the sense of these energies. 

Lemma 32. If S^^{X, J) = M then the identity -Erc(A, J) = Ectm{M) holds. 

Proof. Fill I's, 3's, 5's, . . . into the boxes of the first, second, third, . . . rows of the 
Young diagram A. The sum of these numbers in the a-th column is (X]a=fc"^ ^'^ 
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rrir, 



because the length of the column is {Yla=k''^a.)- By summing them for 1 < a < /i 
we obtain ^^>. -^^ min(i, j)mimj. Then by adding the sum of all the riggings we 
have EYic{)\,J)- Let M = S~^{X,J). Recall the map Uja = ^^ o Qja in (5) and 
consider what it does on the first row of Ma+i{X, J)- The Qjo. embeds J", . . . , J^^ 
somewhere and the $1 adds 1,3,5,... to the first, second, third, . . . columns. By 
summing them for 1 < a < h we obtain Ectm{M) that is certainly equal to the 
above _Erc(A, J). □ 

Let us introduce the notion of partition functions for the box-ball system with a 
specified soliton content [7]. By means of the sets of automaton states (9), (10) and 
the energy in Definition 30 we define 



MeW(L,A) 
MgW+(L,A) 



(15) 
(16) 



On the other hand we define the partition functions for the rigged configurations by 
using (13), (14) and the energy in Definition 31 as 



(\,j)&nc(L,\) 

F+(L,A):= Y 9''^°^^"'^ 

(\,j)enc+{L,\) 

Recall the definition of the g-binomial coefficient [21] 



(17) 
(18) 



— if ^ 777 <^ 77 

nr=i(i-'?'=)nfe;r(i-'?'=) - - ' 

otherwise. 



Then we have that 

Lemma 33. The following relations hold, 

F(L, A) = g-l^l+'^W J] 



i>l 



F+(L,A)=g^Wn 



i>l >- 



Pi{\) +mj(A) +i 
mi{\) 

Pi{\) +mj(A) 
mi{X) 



Here Pi{X),mi{X) and (j){X) are given by (11), (12). 
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Proof. The latter relation follows from the identity of the g-binomial coefficient [21] 

p + m 
m 



,Ji+---+J„ 



E r 

0<Jl<-<Jm<P 

The former relation also follows from this identity by replacing p hj p + i and J/s 

by Jj + i's. 



D 



By Lemmas 28, 29 and 32 the identities Z{L, A) 
Y^{L, A) hold. Thus we have that 

Theorem 34. The following relations hold, 

'Pi{\) +mi{\) +i 



Y{L,\) and Z+{L,\) 



M€n{L,X) 



Ectm(M) 



-|A|+<^(A) 



n 



mi(A) 



Bctm{M) 



>{A) 



n 



'pi{\) +mi{\) 
mi{X) 



M£H+{L,\) i>l 

These are the identities proposed in [7]. 

By taking a sum over all the soliton contents with a fixed number (s = |A|) of 
balls their left-hand sides yield simple expressions (See Appendix A). Then we have 
that 

Proposition 35. The following relations hold. For < s < L 



E'*'n 



Xhs 



i>l 



Pi{\) +mj(A) +i 
mi{\) 



and forO<s< L/2 



L 

s- 1 



E^^TI 



Ahs 



j>l 



mi{\) 



The latter is a g-analogue of the Bethe's formula [22]. In the context of Bethe 
ansatz the summation variable A = (l™i2™2 . . .) has meant that the associated eigen- 
state of the Heisenberg magnet is given by a set of variables which satisfy an algebraic 
equation (the Bethe equation), where the multiplicity of the strings (particular con- 
figurations of its roots in the complex plane) of length i was given by rrij. Thus 
we obtained a new interpretation of a well-known fermionic character formula; the 
summation variable A is representing the soliton content of a cellular automaton. 

On the other hand the former formula may be new and the author expects that 
it can also be used in a completeness problem of Bethe ansatz analysis for some 
quantum mechanical system. 
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9 An alternative description of the inverse scat- 
tering transform 

In this section we describe sl{2) case of the bijection in [8] in our setting and prove 
that it is equivalent to our inverse scattering transform. 

Let (A, J) be the rigged configuration given at the beginning of Section 6. We 
show a procedure to make another rigged configuration (A', J') from it. Let s be 



the smallest integer with rris > that maximizes J^] = J^^ + 2 ^ • ^ min(s, j)mj. 
Call the row of the diagram A of length s and with rigging J^^ the shortest singular 
row [8]. We first delete the right-most box of the shortest singular row, and then 
rearrange the order of the rows of the diagram so that the result is again a Young 
diagram. We replace the rigging of the shortest singular row by a specific manner; 
explicitly 

1. Ifs = llet A'= (l™i-i2™2.../i'"ft), J'l = Ji\{J^J, and J'* = jHoT2<i<h. 

2. If s > 2 let A' = (i"! . . . (s - i)"^s-i+igms^i ...h""^), j"-^ = .r^^ u {.r^i - 1}, 

J'« = .r\ {.r„J, and .P = .r ioii^s,s-l. 

Here J^l = J^^ + 2(m^ H h rrih). 

Theorem 36. If s = 1 the matrix M{\', J') is made out of M(A, J) by removing 
its right-most column. If s > 2 it is made out of M(A, J) by subtracting 1 from its 
lower right element. 

Proof. For both cases we have Ma(A', J') = Ma{X, J) for s + 1 < a < h. Since J'"^ = 
J* \ {Jm^} we see that the matrix Qj's(Ms+i{\', J')) is made out of r2js(Ms+i(A, J)) 
by removing its right-most column '(J^^, Jm,)- By applying $1 we can deduce that 
Ms{\', J') is made out of Ms{X, J) by removing its right-most column *( J^^ — 1, JmJ- 
Hence the s = 1 case follows. Suppose s > 2. Then since J'"^^^ = J"^^^ U { J^^ — 1} we 
see that ^jis-i[Ms{X', J')) and fljs-i{Ms{X, J)) have the same number of columns; 
we also see that the right-most column of the former is XJ^l ~ ^^^ml " l)' ^^^^ 
of the latter is '(J^^ — 1, J^O^ ^^^ ^^^ other columns are equal to each other. By 
applying $1 again we obtain Ms-i{\',J') with its lower right element J^!"^ ~ 1) 
and Ms_i(A, J) with its lower right element J^l~^- By the proof of Lemma 24 one 
can deduce that their difference is kept untouched until we obtain M(A', J') with its 
lower right element J^] — 1, and M(A, J) with its lower right element J^\ D 

From this theorem and its proof we have that 
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Corollary 37. The automaton state ^ ^(A', J') is made out of ^ ^(A, J) hy remov- 
ing its right-most ball. 

This implies that our inverse scattering transform for the box-ball system coin- 
cides with the bijection in [8]. In the s/(2) case the latter method (in our terminology) 
determines the positions of the balls by repeating the procedure for obtaining (A', J') 
from (A, J). We note that although their method was defined only for the highest 
weight states, it can be also defined for the non- highest weight states in the sl{2) 
case. 

Acknowledgements The author thanks Atsuo Kuniba, Masato Okado, and Ya- 
suhiko Yamada for valuable discussions and a collaboration in the previous work 
where he was able to learn the idea of the inverse scattering method. 

Appendix A Partition functions without a speci- 
fied soliton content 

The following items are included here to make this paper to be self-contained, al- 
though they can be found elsewhere. 

Recall the sets of the automaton states (9), (10). We define 

n{L,s) = ux^sn{L,x), 
n-^{L,s) = ux^sn+{L,x). 

(The notation is a bit ambiguous but no confusion should occur.) They are sets 
of automaton states in the finite interval between walls and L, without a speci- 
fied soliton content but with the number of balls set to be s. With the energy in 
Definition 30 we define 

Z{L,s):= J2 g''^™^'^\ 

M<=H{L,s) 

Z+{L,s):= J2 q'"''™'-^^- 

MeH+(L,s) 

Then we have that 

Lemma 38. The partition function Z{L, s) satisfies the recursion relation 

s 

Z{L, s) = Z{L -l,s) + J2 q^''"Z{L -k-l,s-k). (A 1) 

k=l 
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The Z^{L, s) also satisfies the recursion relation of the same form. 

Proof. Consider the array of L boxes between walls and L. The first term in the 
right-hand side of (A 1) is for those automaton states without ball in the right-most 
box. In the summation, the k-th term is for those states which have their right-most 
empty box on the left of the wall L — k. D 

Each of the partition functions is uniquely determined by the recursion relation 
of the form (A 1) under the boundary conditions 



Z{L,0) = Z{L,L) = 1, 
Z+(L,0) = l,Z+{2s-l,s) 

Proposition 39. The following identities hold. 



0. 



Z{L,s) 

z^ms) 

Proof. Using the identities of the g-binomial coefficients [21] 



L 

s 

LI r L " 

s s — 1 



(A 2) 
(A3) 



(A 4) 
(A 5) 



'L 




'L-l 


+ g^~' 


"L-l" 




'L-l 




"L-l" 




= 






= q' 




+ 




s 




s 




s - 1 




s 




s - 1 



one can deduce that the expressions in the right hand sides of (A 4) and (A 5) satisfy 
the recursion relation of the form (A 1). They also satisfy the boundary conditions 
(A 2) and (A 3) respectively. D 

The right hand side of (A 4) is the generating function for the number of partitions 
of an integer into at most s part, and each part is less than or equal to L — s [21]. 
An energy preserving bijection between 7i{L, s) and the set of all such partitions is 
given as follows. Given M G 7i(L, s) consider the automaton state associated to M. 
If its left-most s boxes are filled, then it is mapped to 0. Suppose the state is not the 
case. Then the state admits the following description: there are ai filled boxes from 
the left, then Pi empty boxes, 0:2 filled boxes, P2 empty boxes, ..., Up filled boxes, 
Pp empty boxes. Here p is an integer > 2, X]fc=i '^k = s, and X]fc=i Pk = L — s. We 
assume a^, /9fc > 1 except ai, Pp > 0. Then the state is mapped to the largest Young 
diagram in Figure 3. 
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/9l ^2 



Ps 



pp-1 (3p 



«! + 1 

as 



a«-i 



a,- 



Figure 3: A graphical representation of the energy preserving bijection between 
automaton states to restricted partitions. The largest Young diagram (thick line) is 
made oi p — 1 border strips (skew Young diagrams without 2 by 2 square blocks). 

Conversely given an arbitrary Young diagram we can read off these numbers 
ak, Pk in the following way. We peel off a border strip from the south-east direction 
of the Young diagram. The result is also a Young diagram. Therefore we can repeat 
this procedure as many times as possible. In this way any Young diagram can be 
divided into such border strips uniquely. Then we can read off the numbers ctfc, (3k 
from the lengths of the west and north edges of the border strips. 
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